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The recently proposed Hofava-Lifshitz (HL) theory of gravity is analyzed from the quantum 
cosmology point of view. By employing usual quantum cosmology techniques, we study the quantum 
Friedmann-Lemaitre-Robertson- Walker (FLRW) universe filled with radiation in the context of HL 
gravity. We find that this universe is quantum mechanically nonsingular in two different ways: the 
expectation value of the scale factor (a) (t) never vanishes and, if we abandon the detailed balance 
condition suggested by Hofava, the quantum dynamics of the universe is uniquely determined by 
the initial wave packet and no boundary condition at a = is indeed necessary. 
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I. INTRODUCTION 

In 2009; Hofava proposed a new theory of gravitation fV based on an anisotropic scaling of space x and time t 
coordinates. The resulting theory, since then dubbed Hoi^ava-Lifshitz (HL) gravity, has proved to be power countable 
renormalizable. One of its key points is that, even though it does not exhibit relativistic invariance at short distances. 
General Relativity (GR) is indeed recovered for low-energy limits. Some interesting consequences of this theory 
include the existence of nonsingular bouncing universes and the possibility that it may represent an alternative 
to inflation, since it might solve the flatness and horizon problem and generate scale invariant perturbations for the 
early universe without the need of exponential expansion [SH7] . 

Due to the asymmetry of space and time in the HL gravity, its natural framework is the Arnowitt-Deser-Misner 
(ADM) formalism [8J, where the spacetime metric g^^(t,x.) is decomposed as usual in terms of the 3-dimensional 
metric hij{t,x) of the spatial slices of constant t, the lapse function 7V(t,x), and the shift vector A^*(t, x). In his 
original work, Hofava made an important assumption on the lapse function to simplify the HL gravitational action, 
the so-called "projectablity condition", namely N = N{t). There are, nevertheless, extended models where this 
condition is relaxed. Projectable theories give rise to a unique integrated Hamiltonian constraint, leading to great 
complications when compared with the GR. The so-called non-project able theories, on the other hand, typically give 
rise to a local Hamiltonian constraint, as in GR. Healthy non-projectable extensions of the original HL gravitational 
action are discussed in details in case[S]. Fortunately, since FLRW spacetimes are homogeneous and isotropic, the 
spatial integral can be dropped from the integrated Hamiltonian constraint in our case [ICTl [TT| , yielding a true local 
constraint even for the projectable case. For our purposes here, it suffices to consider the simplest HL projectable 
theory in a FLRW spacetime. We notice that a modified F{R) HL theory in a FLRW spacetime has been recently 
considered in |12| , leading to very interesting results regarding the possible unification between primordial inflation 
and dark energy. Another important assumption originally introduced by Horava is the principle of "detailed balance". 
This condition, which states that the potential in the gravitation action follows from the gradient flow generated by 
a 3-dimensional action, reduces the number of independent coupling constants. Recently, it has became clear that 
the detailed balance condition can be also relaxed |10l [T5HTB] . In particular, in |14l [T5]. the dynamical role and the 
consequences for the matter couplings of the detailed balance condition in classical cosmology are detailed. In this 
paper we will abandon the detailed balance condition, since, as we will show, it gives rise to the most interesting 
quantum universes. The limit where this condition can be recovered will also be discussed. 

There have been many attempts to incorporate Quantum Mechanics into GR. One of the first ones was Quantum 
Cosmology. In Quantum Cosmology, we work with the Hamiltonian (ADM) formulation of GR, use Dirac's algorithm 
|17| of quantization, i.e., the substitution tt^ — >■ —id/6q, where tt^ is the canonical momentum associated with the 
variable q (which can be one of the three canonical variables in GR, hij, N or iV') and the imposition that the 
first-class constraints of the theory should annihilate the wave function of the spacetime. GR has four constraints, 
three of them just tell us that the wave function of the spacetime depends only on the intrinsic geometry of the spatial 
slices in the ADM decomposition, while the last one is a dynamical constraint which gives the dynamical equation of 
quantum cosmology, the so-called Wheeler-DeWitt equation ^18j . The wave function is a priori defined on the space of 
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all 3-metrics. called superspace, which are in general very intricate infinite dimensional spaces. However, we can take 
advantage of the symmetries of a homogeneous universe to freeze out all but a finite number of degrees of freedom of 
the the metric and then quantize the remaining ones. These models are known as minisuperspace models. Quantum 
cosmology in FLRW minisuperspace filled with a perfect fiuid has been shown to be viable and interesting in the 
sense that the initial big-bang singularity is not present in such model since (a) (i) ^ for all times and the classical 
behavior of the universe is recovered for large times |19[ [20]. Moreover, in this class of models, a certain evolution 
parameter of the fluid gives us a measure of time and one can investigate the evolution of the scale factor as the fluid 
evolves. 

For static spacetimes, Horowitz and Marolf [2T] found an original way of classifying a spacetime as quantum 
mechanically nonsingular. In their work, a spacetime is said to be quantum mechanically nonsingular if the evolution 
of quantum particles in the classical background is uniquely determined by the initial wave packet, i.e., no boundary 
conditions at the classical singular points are necessary. This is equivalent to say that the spatial part of the wave 
equation is essentially a self-adjoint operator, i.e., it has a unique self-adjoint extension (for a review about the 
mathematical framework necessary to define quantum singularities in static spacetimes, see (S^)- In the GR context, 
the quantization of the FLRW minisuperspace filled with a perfect fiuid does require a boundary condition at a = in 
order to assure the self-adjointness the Wheeler-DeWitt equation, which, on the other hand, is necessary to guarantee 
a unitary time evolution. Mathematically, the Hamiltonian operator corresponding to the evolution equation of the 
universe is not essentially self-adjoint. In this way the quantum dynamics of the universe is not unique since we do 
not know, in principle, which boundary condition we must apply at the initial singularity. However, as we will see 
later, it is possible to find quantum cosmologies in the HL gravity context for which the quantum evolution of the 
universe is unique, and no boundary condition for the wave function is indeed necessary. 

In this paper we will apply the machinery of quantum cosmology to the HL theory of gravity. In particular, we will 
investigate the necessity of initial boundary conditions for the Wheeler-DeWitt equation and also the behavior of the 
universe, by examining the time evolution of the expectation value of the scale factor. A certain evolution parameter 
of the radiation filling the universe will play the role of the time coordinate. The content of the universe will be 
introduced in the gravitational action via the Schutz formalism [23] [21] , demanding the recovering of the usual GR 
formulation in the low-energy [25j. The paper is organized as follows. Sections II, III, and IV present brief reviews 
of the results we need, respectively, the main definitions about quantum singularities, the HL theory of gravity, and 
the usual quantum cosmology in the GR context. Our main results are presented in the sections V and VI. The last 
section is devoted to some concluding remarks. 



II. QUANTUM SINGULARITIES 

Typical solutions of the Einstein field equations are known to exhibit singularities. They can be classified as |26| : 
quasiregular singularities, where the observer feels no physical quantity diverging, except at the moment when its 
worldline reaches the singularity (for instance, the conical singularity of a cosmic string); scalar curvature singularities, 
where every observer approaching the singularity experience diverging tidal forces (for example, the singularity inside 
a black hole and, more important in the present context, the big bang singularity in FLRW cosmology); non scalar 
singularities, where there are some curves in which the observers experience unbounded tidal forces (for example, 
whimper cosmologies). It is well known [27, that under very reasonable conditions (the energy conditions), which 
basically state that gravity must be attractive, singularities are inevitable in GR. In this way, cosmological models 
with non-exotic fiuids, as radiation or dust, present typically an initial singularity, known as the big bang singularity. 
Since we cannot escape this fact in GR, we hope that the quantum theory of gravitation will solve this issue, guiding 
us on how to deal with the singularities, or even excluding them at all. Unfortunately, we do not have such theory yet. 
However, there are several evidences that this theory would actually solve this problem. These evidences come with 
the introduction of quantum mechanics in GR in many different ways. In this paper we will highlight two distinct 
approaches. 

The first approach is quantum field theory in curved spacetimes. In this framework, we analyze the behavior 
of quantum particles (or fields) in a classical curved background, which we assume to be a regular solution of the 
Einstein field equations. We adopt Horowitz and Marolf's definitions [2T] . In their work, they analyze the behavior 
of a scalar particle in singular static spacetimes possessing a timelike Killing vector field In such spacetimes, the 
wave equation can be separated into 

where A = —VD'^{VDi) + V'^M'^ and V = — ^^C^i, with Di being the spatial covariant derivative in a static slice 
E not containing the singularity. In principle, the domain 2?(A) of the operator A is not known, so we choose as a 
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first attempt X'(A) ~ C^(S). In this way, our operator is symmetric and positive definite. Howerver, this domain is 
unnecessarily small, or in other words, the conditions on the functions are so restrictive that the operator A is not 
self-adjoint. Its adjoint operator A* has a much larger domain 'DiA') = \^ E -^^(S) • G ^^(^)}- It is important 
to notice that we have chosen L^(S) as the Hilbert space of our quantum theory (for a discussion about this point 
see [22l ) . We must relax the conditions on the allowed functions in order to extend the domain of A in such a way 
that T>{A*) — > ^{A). If the extended operator is unique, A is said essentially self adjoint and its extension is given 
by {A,'D{A)), where A is the closure of A (for more detais see [IS]). The time evolution of the particle will be then 
given by 

*(i) =exp(-ia'^')*(0), (2) 

and the spacetime is said to be quantum mechanically nonsingular. However, if the extension is not unique, i.e., 
if there exists infinitely many extensions Aa, with a being a parameter such that to each a there corresponds one 
boundary condition at the singular point, then we have a different time evolution 

«'„(t) =exp(-itAy2)*(0) (3) 

for each a. In this case, the spacetime is said to be quantum mechanically singular. Similarly to the classical case, 
when a spacetime is quantum mechanically singular, an extra information (a boundary condition) must be given in 
order to obtain the time evolution. In GR, in particular, we do need to tell what happen to the particle when it 
reaches the singularity. 

The second approach we exploit here is quantum cosmology in minisuperspaces. In this framework, we consider a 
few degrees of freedom of the system (the rest is assumed to be frozen) and quantize the constraints of the theory via 
Dirac's algorithm. We impose [fljPj,] = i (in units where h = 1), where a is the scale factor of FLRW models, and 
[T,pt] = i, where T is a parameter associated with the evolution of the fluid filling the universe, obtaining, in this 
way, the Wheeler-DeWitt equation of the universe, which, as we will see, is a Schodinger like equation, from where 
we can define an internal product between two solutions and, therefore, evaluate expectation values of observables. 
In this context, we define the universe as nonsingular if (a) (t) ^ for all times. Since the operator d is positive in 
L^(0, oo), we will have (a) (t) = if the wave function representing the universe is sharply peaked at a = 0. Note 
that this criterion is different from that one originally stated by DeWitt, which says that the universe is quantum 
mechanically nonsingular if ^(a — 0,t) ^ V t. In fact, it was shown that this criterion is not enough to prevent 
singularities in quantum cosmological models [29] . 

The two classifications of quantum singularities described above belong to completely different frameworks, but we 
can apply the mathematical machinery used in static spacetimes in order to decide if the evolution of a wave packet 
governed by the Wheeler-Dewitt equation is unique in a given quantum cosmology scenario. 



III. HL GRAVITY 



In order to introduce the HL theory of gravity, let us first introduce the decomposition of the metric in the ADM 
form 

ds^ ^ -N^c^dt^ + h,j{dx' - N'dt){dx^ - N^dt), (4) 

and then let us postulate that the dimensions of space and time are (in units of momentum) [da;*] = — 1 and [dt] = —3. 
This assumption assures that theory is power-countable renormalizable in four dimensions. In these units, we have 
[N] = [hij] — 0, while [A^*] = 2, leading to [ds'^] = —2. Notice that the volume element, defined by 

dVi = NVhd^yidt, (5) 

has dimension [dVi] = —6. 

The extrinsic curvature tensor, which measures how the spatial slices in the ADM decomposition of spacetime 
curves with respect to external observers, is defined by 

It is easy to see that it has dimension [Kij] = 3. The most general term involving the extrinsic curvature tensor which 
is invariant under the group of diffeomorphism of the spatial slices will define the kinetic term in the action. This 
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term depends on two coupling constants a and A and is given by 

SK=aj dt<f^^/hN {K,jK'^ - XK^). (7) 

Note that [a] =0, i.e., a is a dimensionless constant. This is the reason why we made the choice [dt] = —3. 
The potential term for the gravitational action is given by 

Sv ^ - [ dtd^xVhNV[hij], (8) 



where ^[/lij] is built out of the spatial metric and its spatial derivatives. Since [dV4] — —6, we must have [V^[ft.ij]] = 6 
in order to assure that Sy be a scalar. The most general action (without the detailed balance condition) containing 
terms with dimensions less or equal than 6 is given by (for more details see |10| ) 

Shl — Sk + Sv, (9) 

where 

+ 9iR^ + 95R{R^JR"^) + 96R'jR\R'"^ (10) 

Here the constant C has dimension [(,] — 1 and ensures that all the coupling are dimensionless. In order to restore 
the units where c = 1, i.e., [dx] = [dt] we need to perform the transformation dt — !■ C,^^dt. 
Since [Rji^i] = 2. as we go to lower momenta, the dominant action is 

SiH = J dtd^xN^/h [a {K,jK'^ - \K^) - giC^R - goC^] . (11) 
We can now re-scale time and space so that a = 1 and gi = —1, and set c — X — 1, leading to 

Sir = y dtd^xNVh [{K,,K'^ -K^)+R~ 50C'] ■ 



(12) 



Note that, by choosing 



167rG' 

we have the usual Einstein-Hilbert action 



A=^, (13) 



1 

~ 167rG 

where g^^i/ and ^^'>R are the spacetime metric and Ricci scalar, respectively. 
The full HL action we will consider hereafter is 



^ / dtd^TiN\/h (K,,K'^ -K^ + R-2A) 
ttG J ^ ' ' 

*a;^/^(^((4)i?-2A), 



(14) 



. Ml 

'^HL — 
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j dtd^xNVh {K,jK'^ - XK'^ + R-2A 

g2Mp^R^ - gsMp'^R.jR'^ - gAMp'^R^ (^5) 
g^Mp^R{R,,R'^) - geMp^R'^R\R\ 
gjMp^RV^R - gsMp^V.R.kV'R''') , 



where Mp = I/VSttG stands for the Planck mass in c = 1, = 1 units. 
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IV. QUANTUM COSMOLOGY IN GR 

In the so-called Schutz formalism ^51 for the matter content of GR, the four-velocity of a perfect fluid is 
expressed in terms of six potentials in the form 

U, = ^i-\(j),, + a(3^, + 9S^,), (16) 

where /z and S are, respectively, the specific enthalpy and the specific entropy of the fluid. The potentials a and (3 are 
connected with rotations and, hence, they are not present in the FLRW universe due to its symmetry. The potentials 
(j) and 9 have no clear physical meaning. With the usual normalization 

U-'U, = -1, (17) 

Schutz showed that the action for the fluid in GR is given by 

Sf^ f d^x^p, (18) 



where p is the pressure of the fluid, which is related to the density by the equation of state p = wp. In this way, the 
total action for the spacetime flUed with a perfect fluid is given by 

J d^xV^ {R - 2A) + J d^x^p. (19) 



Varying the above action with respect to the metric we get the usual Einstein equations 



G^, + Ag^,^Mp^T^,, (20) 



with r^i, given by 

(21) 

For the FLRW universe with metric 

ds^ = -N^df + a(t)2 ( -f r^dnA , (22) 
\ 1 — fcr^ / 

where dVi^ is the metric in the unit sphere and k = —1,0, 1 for the open, flat and closed universe, respectively, the 
four velocity of the fluid is given by = iV(S° so that 

M = (<^ + eS)/N. (23) 

On the other hand, by thermodynamical considerations, Lapchinskii and Rubakov |30| found that the expression for 
the pressure is given in terms of the potentials by 

^ __UHJ_^_ s/w 
^ (l + w)i+i/» 

/ . . X 1+1/w (24) 



(l+W7)l+l/» I N 

For the particular case of FLRW universes, we have 



a 



2k, 



^ij — h-ij^ Rij — r, hij ^ (25) 

Na 



with hij — diag (yj^kr^'' ''^ ^^"^ ^) ' ^-^^^^ ^'^^'■^ total action is given by (in units where IQttG — 1) 

S = j dtd^xNVh {K,jK'^ -K^ + R) + J dtd^xNy/hp 



d^x / dtl - 6— -I- QkNa 
Vl - fcr2 j 1 N (26) 



6 



The spatial integration does not affect the equations of motion, and we have the following canonical momenta asso- 
ciated, respectively, to the dynamical variables a, (p and S 

N ' P"^^ + ' ' (27) 

PS = dp^,. 

The Hamiltonian of the system will be given by 

H ^Paa + p^{d^ + eS)- L, (28) 

where 



N {\ + wY+^l^ (29) 

X ((/) + 6i5)i+i/'"e-^/'". 

After a tedious but straightforward calculation, we find 

H ^ N {-^-Qka+ pi+™a-3"'e-^/"'^ . (30) 

Since the action does not depend on TV, we conclude that N is actually a Lagrange multiplier of the theory. This 
is not surprising since the results could not depend on how the spacetime is sliced. Varying the action 



S= dt 



Pad + p4^ + 9S)-H (31) 



with respect to N leads to the super-Hamiltonian constraint 



2 

n = -6ka+ pJ,+"'a-3'«g-5/«, ^ q_ .32) 

24a ^ ' 



Performing a canonical transformation of the form 

PS _ ' (33) 



rr S (l+w) s 



(j) ^ (j) + {I + w) 

P(j> 

we get 

2 

24a a^™ ^ ^ 

Now, we proceed with Dirac's algorithm of quantization of constrained systems by making the substitutions pa — >■ 
—id /da, pt — —id/dT, and demand that the constraint annihilate the wave function, finding the Schrodinger- 
Wheeler-DeWitt equation of the universe 

— + 144fca2* + i24ai"3"' — = 0, (35) 
aa-^ ot 

with t = —T being the time coordinate in the gauge N = a^^ , as follows from Hamilton's classical equations of motion 
[31| . Notice that the above equation is of the form id'i'/dt = . In order to the Hamiltonian operator H to be 
self-adjoint we define the internal product of two wave functions as 

/•oo 

= / a^-^'"^*^da (36) 
"'0 

and impose restrictive boundary conditions at a = 0. The simplest ones are the Dirichlet and Neumann conditions: 

*(0,i) = 0, (Dirichlet) 

5*(0,t) ^ , (37) 

= 0, (Neumann). 

aa 

As we will see, the situation is qualitatively different in the HL theory of gravity. 
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V. QUANTUM COSMOLOGY IN HL 

The total action we will consider here is 



S = Shl + J dtd^xNVhp, 



(38) 



where Shl is given by Eq. (15 1. We choose this action basically because the GR action can be recovered in the 



low-energy limit. Discarding the spatial integration again, we have 



5* 



dt 



3(3A - 1 



- 2 

'1^ 



6Nka - 2Aa^ 



12kN 



(3^2+33) 



2AkN 



(954 + 355 + ge) 



(39) 



Let us introduce the following constants (as in Ref. [3^): 

gc = 6fc, 5A = 2A,5r I2k{3g2 + .93), 
5S = 24fc(9(?4 + 3(?5+56)• 



(40) 



Now, proceeding as in GR, i.e., defining the momenta corresponding to each one of the dynamical variables and 
calculating the canonical Hamiltonian, we arrive at 



n 



9 

Pa 



3 , fJr , 9s 

9ca + ffAfl H h ^ 



Pt 

^3w 



0. 



12(3A — l)a a a'^ 

Specializing to the radiation case (w = 1/3), we find the Schrodinger-Wheeler-DeWitt equation 



12(3A - 1) (.gca' - .9Aa' - - ^) * + 12(3A - l)z— = 0, 



(41) 



(42) 



again with t = — T. Note that gr just shifts the energy levels, since it is a constant in the potential. However the term 
gs changes dramatically the effective potential. The case gs — corresponds to the detailed balance condition (see [3] 
where, with the use of the detailed balance condition, Calcagni obtain an action similar to Eq. (39 1, but without the 
term proportional to a~^). 



VI. EXACT SOLUTIONS 
A. Flat FLRW universe 



First, note that if we take a spatially flat universe {k = 0) with A = 0, we have the following equation 

1 _ .d'^ 

"i2(3A- 1) aa2" a^' 



(43) 



which is a Schrodinger-like equation for a free particle with h — 1 and mass mx — 6(3A— 1), except for the requirement 
a > 0. Let us consider only the case A > 1/3. In order to ensure the self-adjointness of the above equation, a boundary 
condition has to be chosen, as discussed in Section III. For sake of simplicity, we choose Dirichlet boundary condition. 
For an initial wave packet of the form 



*(«,0) 



/ 1280-3 



3\ 1/4 



ae 



(44) 



Eq. (43 1 can be easily solved using a specific propagator (see [T5]). The result is 



*(a,i) 



a 



X exp 



m\ + 2it(T ) V 

^a^ \ ( 2itm\a'^c? \ 



-am\ 



m 



A 



40-2^2 J 



exp 



(45) 
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We can calculate the expectation value of the operator a through the formula 

/>oo 








(46) 



Note that if we take A = 1, we recover the result obtained in Ref. ^19j . Nothing changes in HL theory in a flat FLRW 
universe in comparison with GR (as in the classical case, see for instance [5]). In particular, (a) (t) is nonsingular and 
(a) (t) ^ t as t ^ oo, recovering the classical behavior of the universe in the classical limit. Besides, the evolution 
of the wave packet is given once we choose a particular boundary condition at a = 0. Therefore the evolution of the 
universe is not unique in the sense stated in Sec. II. 



B. Closed FLRW universe 



For the spatially closed (fc = 1) FLRW spacetime, there is an effective potential in the Schrodinger-Wheeler-DeWitt 
equation of the universe. From Eq. (42 1, we see that this potential represents a shifted quantum harmonic oscillator 



with a singular perturbation. Setting the mass and frequency of the harmonic oscillator in order to have m\ = 6(3A— 1) 
and uix — ■\/2/(3A — 1), respectively, we have the following Schrodinger-like equation 



1 

2m\ da^ 



dt' 



(47) 



First of all, we will analyze the necessity of a boundary condition at a = on this equation. The first step is to 
separate variables in the form ^'(a,t) = ■tp{a)e^'^^^ , leading to 



Via) 



ip{a) = 2m\Eip{a), 



(48) 



with 



Following Ref. 



V{a) = mluila^ - 2mxgr - {2mxgs)/a^. (49) 
we say that V{a) is in the limit circle case at infinity and at zero if for all A, all solutions of 



+ Via) 



tp{a) = A'0(a) 



(50) 



are square integrable at infinity and at zero, respectively. If V{a) is not in the limit circle case, it is said to be in 
the limit point case. We now enunciate the Theorem X.7 from Ref. [28], which gives us a criterion to decide if the 
Hamiltonian operator H = —d"^ /do? + V(a) is essentially self-adjoint, i.e., if it has a unique self-adjoint extension. 

Theorem 1 Let V{a) he a continuous real-valued function on (0, oo). Then H = —d'^/da^ + V{a) is essentially 
self-ajoint if and only if V{a) is in the limit point case at both zero and infinity. 

We need now a criterion to decide if H is in the limit point or limit circle case at zero and infinity. We will find this 
criterion in the next two theorems, extracted again from Ref. [28J. 

Theorem 2 Let V{a) be a continuous real-valued function on (0, oo) and suppose that there exists a positive differ- 
entiable function M{a) so that 

i) V{a) > -M{a). 

ii) yj M (a)da = oo. 

iii) M{a)/{M{ajf/^is bounded near oo. 



Then V{a) is in the limit point case at oo. 
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Theorem 3 Let V{a) be continuous and positive near a = 0. IfV{a) > near zero then —d? /da?' + V{a) is in the 
limit point case at zero. If for some e > 0, V{a) < (| ^ e)a^^ near zero, then jda? + V{a) is in the limit circle 
case. 

From now on. we will consider < 0. In the end of this section, we will return to the case > 0. Let us first use 
Theorem|2]to show that the Haniiltonian operator in Eq. (47 1 is in the limit point case at infinity. To verify this fact, 



note that the potential ^(a) has a minimum ymin = 

we choose M(a) = 1 and all the requirements of Theorem 



r-^ — 1/2 
I are fulfilled. If Knin < 0, we take Mia) 



at a 



1/4 



,-1/2 



If Knin > 

Klin I- In any 



case, we conclude that H is in the limit point case at infinity. 

Note now that the potential Via) has the form Via) ~ ~2m\gs/a^ near a = 0. Therefore, by Theorem [s] if 
—2mxgs > 3/4 then H is in the limit point case at zero, otherwise it is in the limit circle case. We have established a 
range of parameter in which the operator H is essentially self-adjoint, i.e., if 



-mxgs > 3/8, 



(51) 



then the evolution of the wave function representing the universe is uniquely determined by the initial wave packet 
and no boundary condition at a = is indeed necessary. Otherwise, we need to impose a boundary condition at this 
point. For simplicity we did not consider the case A 7^ 0, but the previous analysis still works in this case. 

It turns out that Eq. (47) can be exactly solved. By introducing the new variable x — y/m\uj\a and a parameter 



Eq. (47 1 becomes 



ipix) ■ 



2 , ^, 
— igr + E) + 
L0\ 4x^ 



i){x) = 0. 



By introducing the new function 2/(77) given by 

V'(x) 



e 2 X 



a+l/2 



it is easy to see that 2/(77) satisfies the associated Laguerre equation 

777/" (77) + (1 + a - 77)?/'(77) + XEviv) = 0, 

where 

'{E + gr) (1 + a)" 



A, 



2a;A 2 

This equation is known to be Hermitian (formally self-adjoint) with the inner product 



if, 9)1 



e '^V°'f{v)giv)dv- 



The general solution of Eq. ( [55| is given by jM| 

1/(77) = AeM{-Xe, 1 + a, 7;) -f BeU{-\e, 1 + a, ??), 



(52) 

(53) 

(54) 
(55) 

(56) 

(57) 
(58) 



where M and U are the confluent hypergeometric functions of first and second kinds, respectively. For \e = n — 
0, 1, 2, . . . , both M and U are polynomials of degree n, proportional to the associated Laguerre polynomial L'^{rj). 
The other linear independent solution is not square-integrable near a = 0, so it must be excluded from the present 
analysis. If ^ N U {0}, we have the following asymptotic behavior for M and t/, as 77 — ^ 00: 



m{-Xe, 1 + «, 77) ^ w-v^e'^r'""""", 

I (-Ab) 
[/(-AB,l + a,77)^7,^^. 



(59) 
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Therefore M is not square-integrable near infinity, meanwhile U is. So M is not an acceptable solution in this case. 
As ?7 — > 0, the asymptotic behavior of U is given by 



U{-\e, 1 + a, 77) - 77-"r(a)/r(-AB), 



(60) 



and, hence, U is square-integrable near a = only if a < 1. For a > 1, we do not have an acceptable solution, 
except in the case A^; = n = 0, 1,2, ... . Therefore, for a > 1 (which corresponds to —rnxgs > 3/8), we quantize 
automatically the energy levels of the universe, which are given by 



En — {2n + 1 + a)uj\ — Qr 
The corresponding normalized eigenstates are given by 



(61) 



r(n + i) 



-,1/2 



/ 2\ 

[mxujxa ) 



_r(a + n+ 1)_ 
-"^a^)L^:(m,^,a^)e-^'^*, 



*„(a,<) =(4toac^a)'/* 
X exp 

and a general solution 5'(a,t), depending on the initial wave packet ^(a, 0), is then given by 

00 

^'(a,i) = ^c„ «-„(&, i), 

with 



(62) 



n=0 



*(a,0)*;(a,0)da. 



(63) 



(64) 



For 1/2 < a < 1, Eq. (60 1 shows that U is square-integrable at a = 0. We then need a boundary condition at a = 
in order to have a well posed Sturm-Liouville problem. They are found in |36| and are given by 



where 



Toy - lim x"+^y\x), 



Tiy = lim 



y{x) + -y\x) 



(65) 



(66) 



li 9 = CO, then lim2,_j.o a:"+^y'(a;) = and we find [37J = n = 0, 1, 2, . . . and the corresponding associated Laguerre 
polynomials. For others values of 9, the quantized energy levels are not so simple. In all cases, the wave function will 
satisfy the DeWitt condition 



*(0,t) = 0, 



(67) 



but we stress the fact that this is not the condition which turns the Wheeler-DeWitt equation into a self-adjoint form. 

It is worth to analyze the case A — >■ 1, — and a — >■ 1/2, where we expect to recover the usual quantum 
cosmology. In this limit, Eq. (62 I becomes 



*n(a,i) 



48 



1/2 



V^22"+i(2n + 1)! 



i/2n+i(Vl2a) 



with 



En = 2n-K3/2. 



(69) 



These same eigenstates have been already found in jT^], in the context of ordinary quantum cosmology. They satisfy 
the Dirichlet boundary condition ^(0,t) = 0. Neumann boundary condition is satisfied when a — ^ ~l/2. We do 



11 



not consider this case here, but it is trivial to generalize our results to — 1 < a < 1/2 (see [SB]). Therefore, the HL 
quantum cosmology tend naturally to usual quantum cosmology in the appropriate limit. 

Have studied the self-adjointness of the evolution equation of the universe, let us now focus on the evolution of the 
expectation value of the scale factor given a solution representing the state of the universe. Obviously, if we calculate 
the expectation value of the scale factor in any of these eigenstates, it will be constant. But the universe evolves so 
that we must consider wave packets representing the state of the universe. In order to find exact solutions, we choose 
an initial wave packet of the form 



*(a,0) = 



2i^+5/2^iy+3/2-| 1/2 



T{iy + 3/2) 



—era 

a e 



(70) 



where gs — ~v{v + \)/{2mx). Note that f/s < in this case, so that the potential V{a) is repulsive, preventing the 
formation of a classical singularity. The propagator for Eq. ( 47 1 is given by [38] 



G(a, a'; t- v) ^'J}^t^2^^^9.t^-i^+^m 
sm [}jj\t) 



X exp 



cot {u:xt){a^ + a'^) 



X J^+l/2 I ^ I 



and through the equation 



\ sin (uJxt) J 



G{a,a';t;i^)^{a',0)da', 



we find, after some tedious calculation. 



2i^+5/2^i/+3/21 1/2 

r(iy + 3/2) 



mxOJx 
i sin {uJxt) 



1^+3/2 



(71) 



(72) 



[2(7 — imxujx cot {LOxt)Y^^^'^ 
X exp 



X exp 



imxtox sin (coxt) cos {10x1)0^ 



2(40-2 sin^ {ujxt) + m\ui\ cos^ {uJxt)) 
(4a2 - mliol) 



The expectation value of the scale factor is given by 



(a) [t) 



2aV {v + 3/2) mxujx 

X \J Aa"^ s\v? {uoxt) + "^a'^a ^^^^ i^xt) 
(a) (0) 



mxuJx 



4o-2 sin {ujxt) + m^ojl cos^ (uJxt)- 



(73) 



(74) 



We stress the nonsingular character of the above equation, since (a) (t) 7^ for all times. As in the usual quantum 
cosmology |19] . the singularity is not present in the quantum model. Note also the similarity with the results found 
in [19_. The only difference here is that the parameter A changes the frequency of oscillation of the scale factor. The 
behavior of this quantum universe in HL theory is not much different from usual quantum cosmology. 

We can also study the case gs > 0, when the potential V{a) is attractive, not preventing the formation of a classical 
singularity. In order to obtain exact solutions we must restrict 2mxgs < 1/8 so that a > 0. We do not have a 
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propagator in this case, but we can use expression (63 1 to find numerically the evolution of the scale factor. This in 
done in Fig. ^ where we have chosen an initial wave packet of the form 



*(a,0) = 2 - 



1/4 



(75) 



The parameters characterizing HL gravity are A = 1, gr = and a = 1/4. Note that (a) (t) 7^ V t. The singularity 
has been excluded. 



A 

O 0.4 
V 



0.2 




FIG. 1: The evolution of the scale factor in the case of an attractive potential. Twenty terms in expansion 1 63 1 were used in 
this numerical approximation. 



VII. CONCLUDING REMARKS 



We have seen that in the HL theory of gravity, it is possible not only to exclude the initial big bang singularity, 
but also to determine uniquely the evolution of the wave function of the universe given an initial wave packet. An 
equivalent statement is that no boundary condition at a = is necessary in a quantum cosmology in the context of HL 
gravity. In general, theories of gravity do not tell us which boundary condition we must choose, so it is a remarkable 
fact that one of these theories excludes this ambiguity. Moreover, in HL quantum cosmology, the evolution of the 
expectation value of the scale factor resembles the evolution found in usual quantum cosmology, the only difference 
being the frequency of oscillation of the bouncing universe. It is interesting to notice that the quantum regime of the 
HL theory of gravity can also provide a viable framework for the description of the "asymptotic darkness" of the visible 
universe |39] . Our early universe results are in a certain sense complementary to the asymptotic regime described in 
|39| . It is certainly worthy to further explore the connections between the two approaches. 
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